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Abstract 

We show that for particles performing Brownian motion in a frozen array of scatterers long-time correlations emerge in 
the mean-square displacement. Defining the velocity autocorrelation function (VACF) via the second time-derivative 
of the mean-square displacement, power-law tails govern the long-time dynamics similar to the case of ballistic mo- 
tion. The physical origin of the persistent memory is due to repeated encounters with the same obstacle which occurs 
naturally in Brownian dynamics without involving other scattering centers. This observation suggests that in this case 
the VACF exhibits these anomalies already at first order in the scattering density. Here we provide an analytic solution 
for the dynamics of a tracer for a dilute planar Lorentz gas and compare our results to computer simulations. Our 
result support the idea that quenched disorder provides a generic mechanism for persistent correlations irrespective of 
the microdynamics of the tracer particle. 
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1. Introduction 

The essence of the Brownian motion of mesoscopic par- 
ticles suspended in a solvent has been understood since 
the pioneering works of Einstein |l]] and von Smolu- 
chowski 12[]. Today, Brownian motion constitutes a ba- 
sic paradigm of stochastic processes ®^ with appli- 
cations in numerous fields, which may even utilize the 
noise as by stochastic resonance |0] and Brownian mo- 
tors 0-01 > not forgetting important generalizations be- 
yond classical statistical physics to quantum Brownian 
system Holl and relativistic Brownian motion Hill . 

Following Einstein and von Smoluchowski, the 
statistics of the particle's trajectories is described in 
terms of a probability distribution for the displacement 
as was measured shortly after by Perrin ll 211 relying 
on the newly developed dark field microscope. These 
displacements are cumulative, i.e., they are the sum 
of many small steps, and provided the steps are un- 
corrected, the central limit theorem applies. Then for 
an unbiased random walk the probability cloud is de- 
scribed by a Gaussian which broadens as time proceeds. 
The corresponding trajectories are continuous but al- 
most nowhere differentiable curves, self-similar in a sta- 
tistical sense. This picture obviously cannot apply to 



the very small time and length scales where collisions 
with individual solvent molecules are resolved. Never- 
theless it seems plausible that on a coarse-grained time 
and length scale all correlations due to these molecu- 
lar processes are rapidly suppressed and Brownian mo- 
tion emerges as mathematically described by a stochas- 
tic differential equation by Langevin 1 1311 with a Wiener 
process as noise term. 

This established picture had to be reconsidered when 
Alder and Wainwright Jl4l[l5ll showed in molecular dy- 
namics simulations that the basic assumption of weakly 
correlated displacements is not fulfilled for the single 
particle motion in a fluid. Rather they found a long-time 
anomaly Br d/2 , B > in the velocity autocorrelation 
function (VACF) due to the conservation of momentum 
in the fluid, where d is the spatial dimension of the sys- 
tem. For a mesoscopic colloidal particle the direct ex- 
perimental observation of these anomalies a la Perrin 
was achieved only 100 years after Einstein's work us- 
ing high-precision photonic force microscopy Jl6Tl . The 
power law correlations persist even in the presence of a 
bounding wall where momentum can be transferred to, 
although they decay more rapidly 1 17- iDll. 



In this paper we argue that quenched disorder consti- 
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tutes a second generic mechanism for persistent corre- 
lations in the velocity autocorrelation function. It was 
shown in the late 60s that in the Lorentz model, where 
a single tracer scatters ballistically in a random array of 
frozen hard obstacles, ring collisions lead to a nonana- 
lytic dependence of the diffusion constant on the den- 
sity of scatterers ll20Tl . The same mechanism implies a 
power-law tail in the VACF JH of the form -Af d+2)/2 , 
A > 0. First, in contrast to Alder's discoveries, the ori- 
gin of this tail is not connected to momentum conser- 
vation, since the tracer exchanges momentum with the 
frozen scatterers. Second, the persistent correlation is 
negative reflecting the caging due to the obstacles. Here 
the correlations in the velocity are inherited from the 
configuration space since, loosely speaking, the particle 
remembers forever which paths are blocked by obsta- 
cles 12211 . The predicted anomalies were soon partially 
confirmed by molecular dynamics simulations 123[ 12411 
yet the prefactor of the tail was by an order of magni- 
tude larger than expected. At higher densities the relax- 
ation appears to become slower suggesting a density- 
dependent exponent |25, 26j]. Later it was suggested 
that a second universal power law takes over as the the 
percolation threshold is approached 12711 in qualitative 
agreement with a theoretical prediction obtained in a 



prean 

mode-coupling approach Il28l 12911 . The scenario of a 
crossover behavior explaining the enhanced prefactor of 
the power law was fully confirmed only recently J30ll . 

Since the kinetic theory developed by Weijland and 
van Leeuwen ll2oTl is rather involved, one is tempted 
to assume that their result constitutes a peculiarity of 
the Lorentz model with few general consequences. In 
this paper we show that the long-time anomalies in the 
VACF persist in a system where particles move diffu- 
sively rather than ballistically through the course of ob- 
stacles. We calculate analytically the first systematic 
correction for the complete scattering function and spe- 
cialize to low wavenumbers to obtain the VACF. The 
long-time tail emerges again due to repeated encoun- 
ters with the same scatterer, yet the derivation of this re- 
sult drastically simplifies compared to the ballistic case, 
since a diffusive tracer finds the same obstacle many 
times without requiring the series of Boltzmann scat- 
tering events as the ballistic particle. We corroborate 
our analytic result by Brownian dynamics simulations 
in the dilute regime and discuss the range of validity 
of the low-density expansion. Last we conclude that 
frozen disorder generically entails long-time algebraic 
decay in the VACF irrespective of the microscopic dy- 
namics implying that the basic assumption of quickly 
decaying correlations in Brownian motion is generally 
not fulfilled in heterogeneous media. 



The notion of universal long-time anomalies also ap- 
plies to hopping transport in disordered lattice models 
A3 lL 13211 where again repeated encounters with the same 
scatterer lead to persistently correlated motion. There 
the VACF has been calculated up to second order in the 
obstacle density and the predictions have been nicely 
confirmed by computer simulation 133 1. 

The paper is organized as follows. In Sec. [2 we 
briefly recall the main results of the anomalies in the 
low-density expansion of the ballistic Lorentz model. 
For the Brownian particle exploring a course of obsta- 
cles introduced in Sec. [3] a multiple scattering expan- 
sion is derived in Sec.|4]revealing that for the first-order 
density expansion it is sufficient to solve the single ob- 
stacle scattering problem. Section [5] calculates the cor- 
responding forward scattering amplitude which is used 
in Sec.|6]to discuss the intermediate scattering function 
and, in particular, the velocity autocorrelation function 
in the dilute case. The analytic predictions are tested 
against computer simulation in Sec. [7] A summary of 
the results and conclusions are given in Sec. [8] 

2. Ballistic Lorentz model and kinetic theory 

In the Lorentz model randomly distributed, possibly 
overlapping obstacles of radius cr represent a random, 
frozen environment in which a single, point-like tracer 
particle moves. The dynamics of the tracer is consid- 
ered as ballistic with elastic scattering whenever an ob- 
stacle is encountered. In particular, energy is conserved 
and the particle's velocity remains constant in magni- 
tude for all times. Diffusion emerges after many col- 
lisions with the scatterer as the direction of the veloc- 
ity is randomized. Since the obstacles are distributed 
randomly and independently, the structures are charac- 
terized solely via the obstacle density n. The only di- 
mensionless control parameter is the reduced obstacle 
density n* := no 4 . 

To lowest order in the obstacle density the motion 
of a ballistic particle in the course of obstacles is de- 
scribed by the Lorentz-Boltzmann equation 13411 . There 
uncorrelated scattering events lead to diffusion on scales 
larger than the mean-free path t ~ cr/n*. In this geo- 
metric problem the time scale is set by the mean col- 
lision time t = €/v by means of the velocity v of 
the particle. The corresponding diffusion coefficient 
Dq ~ v( ~ vcr/n* diverges for small densities since it 
is due to the rare scattering events only that the ballis- 
tic motion becomes diffusive after all. Via a series of 
these uncorrelated scattering events the particle may re- 
turn to an already visited obstacle, giving rise to non- 
trivial correlations. To account for these repeated col- 
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lision events, one has to go systematically beyond the 
Lorentz-Boltzmann equation resulting in a non-analytic 
correction to the diffusion coefficient D. In two dimen- 
sions, d — 2, to which we restrict the discussion in the 
following, one obtains 



Do 
D 



1 



An* 



Inn* - 0.8775n* + 4.519(«*ln«*) 2 + . . 



(1) 

where the Lorentz-Boltzmann diffusion coefficient is 
given by Do = 3v<x/16«*. Uncorrected scatter- 
ing events predict an exponential decay in the veloc- 
ity autocorrelation function Z(f) := (v(f) ■ v(0))/2 = 
(v 2 /2) exp(-4f/3r) with the mean collision rate r _1 = 
2n*v/cr. The repeated scattering from the same obstacle 
then introduces an algebraic tail 



o- 2 1 
Z(f) -8^' 



for t -> oo , n* -> . (2) 



reflecting the infinite memory of the motion in the disor- 
dered system. The amplitude diverges as the system be- 
comes more and more dilute, yet since the initial value 
of the velocity autocorrelation function is density inde- 
pendent, Z(t = 0) = v 2 , the time required to attain the 
tail is set by the growing collision time r. Provided time 
is measured in terms of r, the tail is a small 0(ri*) cor- 
rection to the Lorentz-Boltzmann theory. 

In summary, in the framework of the ballistic Lorentz 
model the non-analytic dependence of D on n* and the 
long-time anomaly emerge in a higher order correction 
beyond the Lorentz-Boltzmann theory. 

3. Brownian particle in a disordered array 

The Lorentz model can be extended to particles that 
perform Brownian motion in the void space, i.e the do- 
main not excluded by the frozen scatterers. In the planar 
case, this might be a useful minimal model for the com- 
plex transport found in cellular membranes or models 
thereof I35H38I1 . The configuration of the environment 
is described by the centers of the obstacles X\,...,Xn 
in a finite hypercubic box of length L. The density of 
the scatterers is kept fixed, n = N/L d , as the thermody- 
namic limit is performed. Due to exclusion, the distance 
of the tracer to any of the obstacles always exceeds the 
radius of the scatterers, |r - x,-| > <x, ; = 1, . . . , N. Pe- 
riodic boundary conditions are assumed throughout for 
convenience. 

A complete statistical description is given in terms of 
the conditional probability density P(r, t\Rt') to find the 
particle at time t at position r provided it was known 
to be at R at an earlier time tf < t. Since the stochastic 



process is stationary, the conditional probability is time- 
translationally invariant and we choose t' — in the 
following. The Smoluchowski equation for P(r, f|R0) is 
then simply the diffusion equation restricted to the void 
space 

d,P(r, f|R0) = D V 2 P(r, f|R0) (3) 

where the nabla operator V acts on the current position 
of the particle r and now Do refers to the short-time dif- 
fusion coefficient. Since the tracer particle cannot pen- 
etrate the obstacles, this equation of motion has to be 
supplemented by the von Neumann boundary condition 

(r - X,) • VP = , for |r - x,| = cr, i=l,...,N 

(4) 

i.e., the flux through the boundary vanishes. 

To avoid difficulties with hard core repulsion, we con- 
sider for the moment a random potential consisting of 
finite spherical barriers, 



U(r) = J] «(r - x,) 



(5) 



where, e.g., u(r) = t/o#(cr - |r|) with the Heaviside step 
function #(■) , and the hard core limit Uq — > oo is antic- 
ipated. Then the propagator obeys the Smoluchowski 
equation 

d t P=-^V-(PVU) + D V 2 P (6) 
k B T 

where k B T is the thermal energy. 

Analytic progress is made by studying the one-sided 
temporal Fourier transform 



G(co; r, R) 



f 

Jo 



e 110 ' P(K|R0)d? 



(7) 



for complex frequencies co in the upper half-plane, 
rather than studying the time dependence of the prop- 
agator P(rf|R0) directly. Then the equation of motion 
translates to 

(-kj - D V 2 )G - -%V • (GVU) = 6(r - R) . (8) 

k B T 

This form of the Smoluchowski equation identifies G as 
an inverse of the Smoluchowski operator and constitutes 
the starting point for the elaborated framework of the 
scattering theory. 



4. Formal scattering theory 

The Smoluchowski equation in the form of Eq. © 
has a mathematical analogy to the time-independent 
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Schrodinger equation, allowing us to employ the tech- 
niques developed for the scattering problem of a quan- 
tum particle. First it is convenient to adopt an ab- 
stract bra-ket notation. The Hilbert space is spanned by 
the (generalized) ket states |r) normalized by (r|r ; ) = 
<5(r - r')- Then we introduce an operator G(lo) with 
a positional representation (r|G(o>)|R) = G(«;r, R). 
The dependence on the complex frequency a> will be 
suppressed in the following. Similarly, we introduce 
the unperturbed Smoluchowski operator Qo via its ma- 
trix elements (r|Qo|^fr) = A)V 2 (r|0) and perturbation 
V by <r|V|*A> = (D /k B T)V ■ ((r\t/r)VU). The Smolu- 
chowski equation © is equivalent to the following op- 
erator equation 

{-ioj - Q)G = I, (9) 

where Q — Qo + V and I is the identity operator. The 
problem is to evaluate the full Green function G in the 
presence of the scatterers. The case of no obstacles 
leads to the unperturbed Green function Go which satis- 
fies 

(-iw-Qo)Go = 1. (10) 

Simple operator algebra reveals that the full Green func- 
tion obeys the Lippmann-Schwinger equation 

G = G + G VG. (11) 

which yields upon iteration the Born series with the first 
order approximation G = Go + GoVGo- The scattering 
matrix defined by S := (I - GoV)~' advances solutions 
of the unperturbed problem to the full one G = SGo- 
It is convenient to single out the event of no scattering 
and to define the T-matrix as GoT = S — I, Then the 
following identities are easily derived 

G = G + G TG , (12) 
T = V + VG T. (13) 

The first relation states that the T-matrix acts like an ef- 
fective potential such that the full multi-scattering pro- 
cess is obtained in first Born approximation. The second 
equation states how the T-matrix is obtained by itera- 
tion, 

T = V + VG V + VGoVGoV + . . . , (14) 

implying that the multiple scattering events are catego- 
rized as single-, double-, . . . scattering events. 

In our case the scattering potential V = v(;) is a 
sum of single obstacle potentials differing merely in the 
position of the scatterer, and thus a multi-scattering ex- 
pansion is appropriate. Let us define the corresponding 
single obstacle t-matrices via 

t© = v(0+v(0G t(0, (15) 



describing multiple scattering form a single obstacle i. 
Then the direct scattering expansion of Eq. ( TBI can be 
reorganized into 

T = Z t; + Z t(/)Got( ^ + Z t (0Got(y)G t(fc) + . . . , 

i i+j i+j 

(16) 

where the first term accounts for the repeated scattering 
with the same scatterer, the second one repeated scatter- 
ing with one scatterer followed by a series of scattering 
events with a different scatterer. The third term con- 
tinues this series including scattering from a third scat- 
terer, which may coincide with the first one. It is now 
clear how the term involving n t-matrices is constructed: 
pick all sequences consisting of n obstacles such that all 
neighboring pairs are distinct. These sequences are rep- 
resented by multiplying the single obstacle t-matrices 
sandwiched with an unperturbed propagator. The mul- 
tiple scattering expansion, Eq. ( TTSI l. remains valid in the 
case of hard obstacles provided the single obstacle t- 
matrix is calculated via the full propagator of the single 
impurity problem analogous to Eq. d!2t . 

The dynamics of the tracer depends on the details of 
its local environment. Performing an ensemble average 
over different initial positions of the tracer and mea- 
suring only the relative displacements reduces the task 
to computing the disorder-averaged propagator G. By 
Eqs. (112113b . this is achieved by equivalently evaluat- 
ing the average many obstacle T-matrix T. The multiple 
scattering expansion shows that to first order it is suffi- 
cient to determine the average single-obstacle t matrix 

T = nL d t + 0(n 2 ) (17) 

since the next terms involve at least two obstacles. 
Furthermore the disorder average restores translational 
symmetry, implying T is diagonal in a plane wave basis. 

The correlations induced by the interaction are con- 
ventionally represented in terms of a self-energy I, de- 
fined via Dyson's equation 

G=G + G IG. (18) 

Comparison with Eqs. (fT2l) . ( TPTI i shows that to first order 
in the obstacle density 

I = nL d t + 0(n 2 ) . (19) 

For the evaluation of the first-order correction in the 
density, it is thus sufficient to solve the dynamics of the 
trace for a single scatterer. 
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5. A single scatterer 

The motion of a single particle diffusing in the presence 
of a fixed spherical hard obstacle is exactly solvable 
which is the main result of this section. Here we use 
an approach which uses a mixed real-space momentum- 
space representation, making the calculation much more 
transparent. Second we choose a formulation reminis- 
cent of the scattering problem of a quantum particle of 
a hard disk using the Hilbert space of the free motion 
as reference. The solution of the quantum scattering 
problem can be found in many textbooks, see for ex- 
ample Ref. |39l Then to account for the exclusion by 
the obstacle a fictitious dynamics of ghost particles is 
introduced, such that the full propagator is a superposi- 
tion of the obstructed motion and a free diffusion pole. 
The hard-sphere scattering also arises in the context of 
the low-density expansion of hard-sphere suspensions, 
solved by Ackerson and Fleishman i|40[] in real space 
and later by Felderhof and Jones 141U43I1 . 

Here we calculate the t-matrix for the single obstacle 
problem via a partial wave decomposition of the solu- 
tion of the Helmholtz equation. To indicate that the op- 
erators evaluated here refer to the single impurity prob- 
lem, we use small letters for the operators. First, it is 
useful to consider the unperturbed Green function go 
(which is of course identical to Go) satisfying 



(-iw - A)V 2 )<r|g |R> = <J(r-R), 



(20) 



without restrictions on the values of the initial and ter- 
minal positions r and R. The propagator g for the 
Brownian particle moving in a two-dimensional plane 
in the presence of a single hard core obstacle of radius 
cr placed at the origin satisfies the equation of motion 

(-ioj - D V 2 Xr|g|R> = 6(r - R)0(|R| - cr) . (21) 

Here the Heaviside step function #(|R| - cr) accounts for 
the fact that the tracer cannot start within the obstacle 
implying that the matrix elements vanish (r|g|R) — 
for |R| < cr. The hard-core repulsion imposes a no-flux 
boundary condition on the propagator, 



r ■ V<r|g|R) = for 



(22) 



at the surface of the obstacle. To make contact with 
the formal scattering theory, the propagator g is to be 
interpreted as the resolvent of some time-evolution op- 
erator. In the case of a hard disk, the dynamics has been 
defined only for particles initially outside of the obsta- 
cle. Completing the description, we assume that parti- 
cles positioned inside an obstacle at the beginning be- 
have as ghost particles diffusing freely without feeling 



the exclusion. The resolvent associated to this operator 
then consists of the propagator we seek plus a diffusive 
pole with a residue corresponding to the packing frac- 
tion of the single obstacle on the plane. 

The solution will be obtained by switching to a mixed 
coordinate-momentum representation 



= Jd 2 R<r|g|R>e i * R , (23) 



<r|g|q) 



where the plane wave states |q) are defined by their over- 
lap with the positional representation (r|q) = e iq r / VA. 
The wavenumbers q form a discrete set such that the po- 
sitional representation obeys periodic boundary condi- 
tions in the plane A — L 2 ; furthermore they are properly 
normalized (k|q) = 5k,q- Then the two Green functions 
satisfy 



{-ia) - A)V 2 )<r|g |q) 



1 : 



tq-r 



VI 



2\/„i„i„\ _ ^ „iqr 



(-k>-D V z )<r|g|q> = — e 
VA 



(24) 
(25) 



where for (r|g|q) satisfies the no-flux boundary condi- 
tion, Eq. (l22l . for |r] > cr and describes the ghost parti- 
cle for |r| < cr. The free motion allows for a plane wave 
solution 

e^V Va 



<r|golq> 



(26) 



-ia> + Doq 2 ' 

which also is the solution for (r|g|q) for terminal posi- 
tion inside the obstacle. From the observation that the 
difference of Eqs. (l24l i and d25l l. 



(iw + D V 2 )<r|g-g |q> = 0, 



(27) 



vanishes for terminal positions outside the obstacle, the 
remaining case |r| > cr follows immediately. The gen- 
eral solution decaying rapidly at infinity of the two- 
dimensional source-free Helmholtz equation ( f2Tb in po- 
lar coordinates is given in terms of modified Bessel 
functions of the second kind K v (-) as 



<r|g - g |q> = J] a m K m ( M ry m * .. 



(28) 



where we abbreviated // = yj-ico/DQ — ( 1 - i) yfaTjlDo 
and ip = Z(r, q). The coefficients Qui £LT6 determined by 
the boundary condition r ■ V(r|g|q) = 0. Noting that 
(r|go|q) can be represented using the Jacobi-Anger ex- 
pansion H Eq. 9.1.41] 



£ i"/„(z)e" 



(29) 
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one obtains 



-1/VA , m qJ' m {qcf) 



(30) 



The propagator can now be evaluated in the momen- 
tum representation via Fourier transform. Here we shall 
need only the forward scattering matrix element 



<ql9 



■golq> = -^= f dW^rlg-g 

yA J\r\>tr 

I 



ilq> 



-2n/A yi qJ'Jqa-) 
ico + D q 2 ^ fiK^cr) 



f 

•J IT 



(31) 



rdr J m (qr)K m (fir) , 



where we used the fact that for terminal position in- 
side the obstacle, the free and full propagator are identi- 
cal. The second line is obtained using again the Jacobi- 
Anger expansion and integrating over the relative angle. 
With the help of the indefinite integral 



rdrJ m (qr)K m (jur) 



— r [qK m Qxr)J' m (qr) - /j.K' m (jjr)J m (qr)] , (32) 

q L + /r 



and the relation [44, Eq. 9.1.76]) 

oo oo 

2 J m (zf = 1 => J] JmWjz) = 0, (33) 

m=-oo m=-oo 

one obtains the forward scattering amplitude in closed 
form 



<qlg - golq> = 

2ko- 2 /A 



(-i(0 + D q 2 } 



D q 2 J j; n (qa) 



2 K m (jitr) 
Ho-K'Jjicr) 



(34) 



The t-matrix for the current case is again defined via the 
operator relation g = go + gotgo, and since the unper- 
turbed propagator is diagonal in the wave number rep- 
resentation 



<kl9olq> 



1 



(35) 



-iaj + Doq 2 

the t-matrix for forward scattering is readily obtained 

/ Itl \ 27r °" 2 r, -> V Hi \2 K mQxr) 



Note that the forward scattering amplitude is indepen- 
dent of the position of the scatterer. The self-energy I 
is to first order in the density, n* — cr 2 N/A, equal to 
the average T-matrix of the single scattering t-matrices, 
Eq. ©, 

<q|I|q> = -2mi*D q 2 f j' in(qcT ) 2 ^0ZL + (n* 2 ) 



(37) 



which includes the motion of the ghost particles with 
residue n*n. Since averaging over the disorder restores 
translational symmetry the self-energy is diagonal in the 
wave number representation. The average propagator 
then reads 



<q|G|q> = 



{-ico + D q 2 ) 



27m*D q 2 ^ J' m (qcry 



K m (jxcf) 



IMrK' m (jMf) 



(38) 



up to order 0(n*). This form still contains the diffusive 
motion of a the ghost particle. Yet subtracting the cor- 
responding pole yields up to first order in the packing 
fraction only (1 - n*n) as common prefactor, which re- 
flects the total probability for a randomly placed tracer 
not overlap with an obstacle. Dropping this factor again, 
one can take the preceding result as the conditional 
propagator for particles initially in the void space. 

For reference we give also the corresponding result in 
three dimensions, rC, , = Ncr 3 /I?: 



<q|G|q>=[(-iw + D o9 2 ) 



Ij'Aqc)] 2 



iacrk' f (ficr) 



(39) 



6. Velocity autocorrelation function 



The simplest quantity characterizing deviations from 
simple diffusion is the averaged mean-square displace- 
ment 6r 2 (t) or the velocity autocorrelation function de- 
fined via 



z(f) = idW 



t>0. 



(40) 



The corresponding one-sided Fourier transform Z(w) = 

r-OO 

J e Z(f)d? is related to the long-wavelength limit of 
the Green function G(q,co) := <q|G|q), see e.g. 114511 . 



For completeness the derivation is repeated here. Since 



6 



the mean-square displacement is the second moment of 
the averaged real-space propagator, the VACF can be 
represented as 



The averaged intermediate scattering function 

F(q,t):= j e iqr P(rf|00)d''r , 

exhibits a long wavelength expansion F(q, t) 
q 2 6r 2 (t)/2d + 0(q 4 ), implying 

Z(f) = - lira \d*F(q, t) . 
?-»0 q- 



(41) 

(42) 
1 - 

(43) 



Since G{q,oj) corresponds to the one-sided Fourier 
transform of F{q, f), this last relation can be translated 
to the frequency domain as 



Re Z(oS) = lim — Re G(q, u) , 

9^0 q l 

and to the self-energy Y,(q, u>) = (q|I|q) via 

Z(a>) -Dq - lim — Y.(q, to) . 

q^o q 1 



(44) 



(45) 



Specializing again to the planar case, Eq. ( 1381 ), only the 
terms m = ±1 contribute and we obtain to lowest order 
in the density 



Z(lo) — Do + 7r«*D( 



jicrK' (jicr) 



(46) 



where the frequency dependence is encoded in /u — 
(1 - z) Vw/2Do. For reference let us provide also the 
corresponding result for the three-dimensional case 



it \ r> ±_ An * r> 



(47) 



The velocity autocorrelation inherits the dynamic cor- 
relations due to the frozen landscape of obstacles, in 
particular, it displays non-analytic behavior for low- 
frequencies 



itocr 1 I A-Dq \ \do~j 

1 In r- + 

2D \-iwo- 2 / Do 

+ 0{u> 2 In 2 oj) , (48) 



where y = 0.577 . . . denotes the Euler-Mascheroni con- 
stant. By a Green-Kubo relation, the zero-frequency 



limit yields the long-time diffusion constant D in the 
presence of obstacles, 

D = Z(co = 0) = D (l - nn ) + 0(n 2 ), (49) 

highlighting the suppression of transport due to the ex- 
cluded volume. A Fourier backtransform 



Z(t) 



2 r 

n Jo 



[Re Z(ai)] cos(Lot)da> , 



(50) 



shows that the non-analytic low-frequency expansion 
corresponds to a long-time anomaly, 



z(0 



2fi 



as 



(51) 



with the same power-law as for the ballistic planar 
Lorentz model. Here the prefactor is first order in 
the density reflecting the fact that the tracer can en- 
counter the same obstacle by diffusion many times with- 
out going through a series of scattering events with other 
frozen obstacles. 

The diffusive motion at high frequencies is singular 
too, 



Z(oS) = D 



Dq 



(52) 



which is reminiscent of the skin effect for electromag- 
netic waves on a metal. The corresponding skin pene- 
tration length 6 :- y/2Do / at is the characteristic length 
scale for diffusive transport, and the reduction reflects 
the probability to encounter an obstacle within that 
length. A Fourier back-transform yields the singular 
short-time behavior 



ttD 



Z(t) = -n* A) A — ^ , r^O, (53) 

V <T z t 

characteristic for Brownian motion in an environment 
of hard obstacles. 

7. Simulation results and Discussion 

We have performed Brownian dynamics simulation 
which test our analytic result for the VACF and explore 
the range of validity of the first order approximation to 
the low-density approximation. 

The array of immobilized obstacles is generated by 
placing randomly hard disks of radius cr in a plane of 
size A — Lr . The positions of the scatterers are inde- 
pendently drawn from a uniform distribution and scat- 
terers may overlap which in principle occurs at any den- 
sity n = N/A. To minimize finite-size effects periodic 
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Figure 1 : Dimensionless negative VACF -Z(f) for the dilute planar 
Lorentz gas for Brownian dynamics. Symbols correspond to simula- 
tion results, full lines to the first-order density approximation. Density 
increases from bottom to top. Inset: Rescaling of -Z(f) over n* 



Figure 2: Negative VACF ~Z(t) rectified by the short-time prediction 
(f/fo)" 1 ' 2 . Symbols correspond to our simulation results, full lines to 
the first-order-density approximation. Density increases from bottom 
to top. 



boundary conditions are employed, with typical system 
sizes of L/cr = 10 4 . 

A single tracer explores the void space of the frozen 
landscape by Brownian motion. Special care has to be 
taken to account for the hard-core exclusion at the ob- 
stacles. Here we have relied on an event-driven algo- 
rithm developed for hard sphere liquids ll46ll recently 
employed also for the three-dimensional Lorentz model 
close to the percolation transition 14711 . The basic idea is 
to compute ballistic trajectories including the collisions 
with obstacles, which are interrupted by fictitious colli- 
sions with a solvent acting as a heat bath. In its simplest 
version, these kicks from the solvent are instantaneous 
at regular time intervals with period Tg where a new 
velocity is drawn at random from a two-dimensional 
Maxwell distribution of variance v 2 . At time scales 
large with respect to the algorithmic time Tg, the mo- 
tion is Brownian and solves the diffusion equation with 
a short-time diffusion constant Do = v 2 Tg/4. By con- 
struction, the tracer never leaves the void space and the 
hard-core repulsion is manifested merely in the usual 
specular scattering at the surface of the obstacles. 

The simulations for Brownian tracer particles in the 
low-density range include 5 trajectories for each of the 
155 different obstacle realizations drawn for each den- 
sity, except for n* - 0.01 where 500 different obstacle 
realizations where examined. We measure time in terms 
of microscopic scale to :- cr 2 /Do, i.e., the time needed 
for the particle to diffuse one obstacle radius without 
obstruction. The algorithmic time tb should be much 
smaller than to and here we used r# = 0.0025fo- The 
negative velocity autocorrelation function -Z(f) is dis- 



played in Fig.[T]on double-logarithmic scales covering 
four non-trivial decades in time and more than three or- 
ders of magnitude in signal. For the smallest density 
n* = 0.01 the data coincide with the theoretical first- 
order density approximation for all times. The time axis 
in the figure starts at t = Tb and the small increase in the 
first data points visible is still affected by the algorithmic 
resolution. The curves corresponding to moderate den- 
sities still follow the first-order low-density prediction at 
short times but start to deviate at long times. The long- 
time decay is slower than the expected one quantified 
by an apparent exponent which becomes smaller upon 
increasing the obstacle density. The crossover regime 
shows only a slight flattening of the curvature consistent 
with the observed increase of the exponent. The inset in 
Fig.[TJcorroborates very nicely the direct density depen- 
dence of the VACF following from Eq. (1481) . All simu- 
lation results overlap in the rectification with the theo- 
retical curve besides the exponent variations discussed 
above. Let us mention that for the ballistic case a nu- 
merical confirmation of the long-time tail with univer- 
sal exponent has been achieved only recently ifjoll . yet 
the amplitude for the power law still deviated by 25% 
from the theoretical value for n* = 0.005. Furthermore 
our simulation for a Brownian walker allows to compare 
the VACF for all times to the theory since the VACF is 
known beyond its asymptotic behavior. 

For a more thorough examination of the short-time 
behavior (Fig. |2]i it is advantageous to rectify the VACF 
with respect to the predicted f -1 ' 2 behavior of Eq. d53V 
For the smallest density again a perfect agreement is ob- 
served. The moderate densities are still well described 
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Figure 3: Long-time behavior of Z(i) rectified by the long-time tail 
(t/to)~ 2 . Symbols correspond to the simulation results, full lines are 
the first-order-density approximation. Density increases from bottom 
to top. 
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Figure 4: Long-time diffusion coefficient D as a function of reduced 
obstacle density n* = no 2 . Symbols correspond the simulation re- 
sults, the straight line D/Dq = 1 - nri" is the first-order-density ap- 
proximation. 



by the first-order-density theory reflecting the fact that 
the particle did not have time to undergo collisions with 
more than one obstacle. At short times the effects aris- 
ing from a finite tb amplify as the system becomes 
denser, because the number of collisions per Tg is aug- 
mented by up to a factor of 15. For the highest density 
n = 0.15 we measure already 0.026 collisions on aver- 
age between the algorithmic assignment of new random 
velocities. Assuming a Poissonian distribution for the 
scattering events the probability that two collisions take 
place in this time interval is roughly about 0.1% and 
therefore not completely negligible. 

A rectification plot for the long-time behavior is de- 
picted in Fig. [3] The VACF reaches the predicted scal- 
ing r 2 for all densities, yet the observed plateau val- 
ues at long times increase stronger than expected. For 
the ballistic Lorentz gas it has been shown that such 
an increase arises due to a competition of the long- 
time tails and the critical relaxation at the percolation 
transition [f30h - There, the underlying fractal structure 
induces anomalous transport for the mean-square dis- 
placement [48] resulting in a fractal power-law decay of 
the VACF. The dynamics of the two-dimensional Brow- 
nian motion close to the percolation transition shall be 
discussed elsewhere l49ll . 

The long-time diffusion coefficients were extracted 
via D - lim,^ 00 (l/4)d(5r 2 (f)/df from the simulated 
mean-square displacements. The suppression of the 
long-time diffusion coefficient predicted in Eq. d49l is 
shown in Fig. |4] Up to densities of about n* = 0.2 
the data follow quite nicely the predictions; above, 
higher order corrections gain importance. Neverthe- 
less the first-order theory provides a remarkably good 



description of the data for all densities. Extrapolation 
suggests that the diffusion coefficient should vanish at 
n* = 1/tt = 0.318.. which is surprisingly close to the 
measured critical density of the percolation transition 
n* = 0.359.. IU. 

8. Conclusion and Outlook 

The notion that correlations quickly die out at time 
scales beyond some characteristic relaxation time of the 
problem has been shown to be incorrect in general. Be- 
sides the meanwhile established long-time anomalies 
due to momentum conservation in fluids Jl4l[l5Tl . a sec- 
ond paradigm leading to persistent correlations is iden- 
tified. Quenched disorder implies repeated encounters 
with the same obstacle, and the information encoded in 
the exclusion of the configuration space manifests itself 
in measurable quantities such as the mean-square dis- 
placement or the velocity autocorrelation function. 

Previous studies focused on the ballistic motion in 
quenched disorder Ir20[ 21] where the theoretical analy- 



sis is quite involved, and the long-time anomalies could 
be considered as merely a peculiarity of the Lorentz 
model. The identification of similar persistent correla- 
tions for hopping transport in disordered lattices i3ll - 
3311 suggests that memory effects may apply to a much 
larger class of systems. Here we have calculated the 
memory effects for a Brownian particle in a random en- 
vironment of hard scatterers to first order in the obstacle 
density. We conclude that the only ingredient necessary 
for the long-time tails is the frozen disorder. Since dis- 
order is ubiquitous in nature and the effects arise at all 
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obstacle densities, we conclude that the power laws in 
the VACF are present for all real systems. 

We have confirmed our analytical results by computer 
simulation for Brownian motion in a disordered array 
of obstacles. The Brownian particle follows the first- 
order theory quantitatively at low densities and qualita- 
tively for moderate ones confirming that the long-time 
anomaly persists for all densities where diffusion occurs 
in the long-time limit. Furthermore, we have shown that 
the first-order-density expansion gives a reliable picture 
for Brownian dynamics, in contrast to the ballistic case. 

The Brownian dynamics in the presence of hard ob- 
stacles displays a second power law in the VACF at short 
times which is due to single scattering events from an 
obstacle. These effects have no analog in the ballistic 
case nor in hopping transport in disordered lattices. 

In the present study the unobstructed dynamics is 
considered to be a random walk where momentum does 
not play a role. For a colloidal particle suspended in a 
fluid and moving through a dilute course of obstacles, 
one may expect that both the vortex diffusion of mo- 
mentum in the fluid as well as the repeated encounters 
with static obstacles leads to an algebraic long-time de- 
cay of the VACF. Since the decay t~ d ' 2 due to vortex dif- 
fusion for unconfined motion is slower than the one for 
the obstructed motion r (d+2) ^ 2 one may anticipate that 
momentum conservation is more important than steric 
hindrance at least for long times. Yet, the obstacles can 
also carry away momentum and the vortex diffusion in 
the fluid should be cut off at time scales where the agi- 
tated fluid encounters an obstacle. For the case of a wall 
in three dimensions it has been shown that the long- 
time anomalies are suppressed 117-3 to r 5 ^ 2 rather 
than r 3 ^ 2 . Hence, for a disordered system in a fluid a 
competition between both scenarios appears to be rele- 
vant, yet no theory is available. Computer simulations 
for dense hard sphere liquids have indicated a crossover 
from a positive tail due to vortex diffusion as well as 
a negative tail in the densely packed regime, where the 
cages act as a quasi-frozen disordered environment [ 50] . 
This phenomenology has been recently corroborated for 
Lennard- Jones particles J5lll . 

Hydrodynamic fluctuations in two-dimensional fluids 
lead to a series of divergences in the correlation func- 
tions. The diffusion coefficient in an infinite system is 
expected to be (logarithmically) divergent since the fluid 
cannot carry momentum away fast enough. For a dilute 
density of static obstacles one may hope that some of 
these effects are regularized in a natural way, giving an 
intriguing interplay of non-analytic behaviors; a theory 
for such a scenario, however, remains an open question. 
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